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We study the nature of spin excitations of individual transition metal atoms (Ti, V, Cr, Mn, Fe, Co 
and Ni) deposited on a Cu2N/Cu(100) surface using both spin-polarized density functional theory 
(DFT) and exact diagonalization of an Anderson model derived from DFT. We use DFT to compare 
the structural, electronic and magnetic properties of different transition metal adatoms on the 
surface. We find that the average occupation of the transition metal d shell, main contributor to the 
magnetic moment, is not quantized, in contrast with the quantized spin in the model Hamiltonians 
that successfully describe spin excitations in this system. In order to reconcile these two pictures, we 
build a multi-orbital Anderson Hamiltonian for the d shell of the transition metal hybridized with 
the p orbitals of the adjacent Nitrogen atoms, by means of maximally localized Wannier function 
representation of the DFT Hamiltonian. The exact solutions of this model have quantized total spin, 
without quantized charge at the d shell. We propose that the quantized spin of the models actually 
belongs to many-body states with two different charge configurations in the d shell, hybridized 
with the p orbital of the adjacent Nitrogen atoms. This scenario implies that the measured spin 
excitations are not fully localized at the transition metal. 

PACS numbers: 73.22.-f,73.22.Dj 


I. INTRODUCTION 

The Cu(lOO) surface coated with a CU 2 N monolayer 
has turned out to be a remarkable systemic— to probe 
and engineer the electronic properties of individual tran¬ 
sition metal atoms using scanning tunneling microscopy 
(STM) and inelastic electron tunneling spectroscopy 
(lETS). A variety of breakthroughs have been reported 
on this system, such as the first measurement of the mag¬ 
netic anisotropy of an individual quantized spin by means 
of lETSf^ the demonstration of single atom spin torque?^ 
the fabrication of nano engineered chains both with 
antiferromagnetioi and ferromagnetio^i broken symme¬ 
try ground states, probed by means of spin-polarized 
STM, the measurement of spin excitations in spin chains 
with strong quantum fluctuations that prevent spin sym¬ 
metry breakingfii^ the measurement of single spin re¬ 
laxation time by means of voltage pulse pump-probe 
technique)^ the observation of renormalization of mag¬ 
netic anisotropy due to Kondo exchange interactions^ 
and the imaging of spin wave modes with atomic scale 
resolutions^ 

The system has been studied from the theoretical 
standpoint, using a variety of approaches iSiiSrSi Im¬ 
portantly, both the spin excitation spectra of individ¬ 
ual atomsSiSiii and multi-atom structuresas 
well as their spin relaxation dynamics have been suc¬ 
cessfully described using model HamiltoniansiSr— where 
quantized spins interact with each other via Heisen¬ 
berg couplingand are Kondo coupled both to 
the tunneling electron o^^d^ and to the substrate.^SiSSr^i 
Treating Kondo coupling up to second order in per¬ 
turbation theory accounts for spin relaxatio n^^i^^ and 


magnetic anisotropy renormalizationi^i^ Furthermore, 
calculation a^^i^^ up to third order are also able to account 
for non-trivial lETS line shapes, including Kondo peaks. 
Numerical renormalization group non-perturbative cal¬ 
culations for the anisotropic spin Kondo model also pro¬ 
vide very good description for both the finite energy 
spin excitations and the zero bias Kondo peak in these 
systems i^^d^d^ . The origin of the Kondo couplings in 
these systems can be traced down to a multi-orbital An¬ 
derson model for these S' > 1/2 systems^l in line with 
the very well known results for the mapping^^— of the 
single orbital Anderson model^ to the Kondo model. 

In spite of the success of quantized spin Hamiltoni¬ 
ans to describe many experiments of transition metals 
on CU 2 N, there is a problem of principle that we address 
in this paper. Density functional theory (DFT) calcu¬ 
lations, most of them^ i^^i^^~ — dealing with Ti, Mn, Fe 
and Co adatoms on Cu2N/Cu(100), show that nor the 
charge, neither the magnetic moment of these magnetic 
atoms are quantized. 

Here we provide a comprehensive and comparative 
study of the electronic and structural properties of the en¬ 
tire series of 3c? transition metals (Ti, V, Cr, Mn, Fe, Co 
and Ni). To the best of our knowledge, no DFT calcula¬ 
tions have been reported for V and Ni on Cu2N/Cu(100). 
Our calculations conhrm the fractional nature of the av¬ 
erage occupation of the d levels in these systems, which 
is not surprising given their conducting nature, but it 
poses an apparent contradiction with the quantized spin 
model description. In the second part of the manuscript 
we provide a solution to this apparent conflict. We build 
a multi-orbital Anderson model, using as starting point 
the representation of the DFT Kohn-Sham Hamiltonian 
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in a basis of maximally localized Wannier functions. The 
Anderson model includes spin-orbit interactions, crys¬ 
tal field interactions, on-site Coulomb repulsion and hy¬ 
bridization of the d shell of the transition metal with its 
Nitrogen neighbors. We solve the model by exact numer¬ 
ical diagonalization within a restricted Hilbert space that 
includes both d” configurations with configura¬ 

tions, where n stands for the number of electrons in the d 
shell and m stands for the number of electrons in the Ni¬ 
trogen p orbitals. Our numerics show that the low energy 
excitations of the model can be mapped into quantized 
spin Hamiltonians. Within this picture, it is apparent 
that this quantized spin S describes the quantum num¬ 
ber of many-body wave functions that mix states with 
n and n -I- 1 electrons in the d shell, and have thereby 
non-integer average occupation. 

The paper is organized as follows. In Sec. [Ill we de¬ 
scribe DFT calculations for different TM atoms at CU 2 N 
paying particular attention to the structural properties 
(|HBI) . electronic properties (jll Cl) and magnetic proper¬ 
ties (IHDI) . In Sec. IHII we build the multi-orbital Ander¬ 
son model, using as starting point the DFT calculations 
and we analyze the connection with the Spin Models. Fi¬ 
nally, Sec. IIVI contains a summary and a discussion our 
most important findings. 


II. DENSITY FUNCTIONAL CALCULATIONS 
A. Methods 

Most of our DFT calculations of 3d transition metal 
adatoms adsorbed on Cu2N/Cu(100) were done using the 
generalized-gradient approximation (GGA) and GGA-I-U 
for exchange-correlation energy;^ using plane-wave basis 
sets and Projector Augmented-Wave (PAW)^ as imple¬ 
mented in QUANTUM ESPRESSO (QE) code,^ Addi¬ 
tionally, in some particular cases we have performed com¬ 
plementary calculations using Local Spin Density Ap¬ 
proximation (LSDA)^ for exchange-correlation energy, 
using all-electron full-potential linearised augmented- 
plane wave (FP-LAPW) as implemented in ELKi^ 

In order to test the size convergence, we have used two 
super cells with different sizes. Both cells have 4 slabs 
of Cu(IOO), separated by a vacuum region of 15 A. The 
smaller cell has 37 atoms (1 TM, 4 N and 32 Gu ) and 
a bigger supercell, has 82 atoms (1 TM, 9 N and 72 Gu 
). The corresponding structures are shown in Fig. [TJ In 
the smallest structure, the intercell distance between TM 
atoms along the N direction is 7.2A, while in the bigger 
one this distance is 10.SA. 

QE calculations are done in two stages, structural re¬ 
laxation and electronic structure calculation. In the re¬ 
laxation stage, the mesh in k space for the small and big 
cells were 6x6x1 and 4x4x1 respectively. In both cases 
the relaxation was performed until the forces acting on 
atoms were smaller than 10“^ a. m.. In the second stage, 
the meshes used were 8 x 8 x 1 for the small supercell and 



FIG. 1. (Color online) Structure of the TM@Cu 2 N. (a) Small 
supercell where a, the intercell distance between TM atoms 
along the N direction, is 7.2A. (b) Big supercell, with a = 

10.8 A. 


6x6x1 for the big one. In all the calculations we used a 
smearing with a broadening parameter of 0.01 — 0.02 Ry, 
in line with previous workj ^^i^^A^ and we fixed the cut¬ 
off energies for the wave function and charge density at 
40-80 Ry and 400-800 Ry respectively. For the Elk calcu¬ 
lations, we started from the relaxed structures obtained 
with QE. In this case the mesh in k space was 4x4x1, the 
product of the muffin-tin radius and the momentum cut¬ 
off is i?MT fcmax = 6 and we employed Local Spin Density 
Approximation (LSDA)^ for exchange-correlation and 
DFT-I-U with Yukawa screening.— 

In the case of Fe and Go we have also obtained the so 
so called maximally localized Wannier functions (MLWF) 
associated to the Bloch states of the DFT calcula¬ 
tion, using the package WannierQO. The Wannier func¬ 
tions form an orthogonal and complete basis set that we 
can use to describe our system. Importantly, the repre¬ 
sentation of the Kohn-Sham Hamiltonian in the DFT ba¬ 
sis provides an effective tight-binding model to describe 
the electronic states of the system, that we use as a start¬ 
ing point to build a multi-orbital Anderson model, as 
described in Sec. IBIl 


B. Structural properties 

We now discuss the structural properties of a single 3d 
TM atom bonded to the Gu site of the Gu2N/Cu(I00) 
surface. This is the binding site most frequently reported 
in the literature In order to refer to the dif¬ 

ferent TM atoms we shall use indistinctly their chemical 
formula (Ti, V, Cr, Mn, Fe, Go, Ni) or the nominal charge 
on the d shell qd = 2, 3,4, 5, 6, 7, 8. 

CU 2 N is knowii^i^ to form a weakly buckled square 
lattice on top of the Gu(lOO), consistent with our DFT 
calculations. In all cases considered, the adsorption of 
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Surface Cu 


FIG. 2. (Color online) Schematic diagram of the CU 2 N ( Cu 
atoms in brown, N atoms in blue) surface with the TM atom 
(green), dw is the distance from the TM atom to the nearest 
neighbors N atoms, dcu refers to the distance from the TM 
atom to the Cu atom lying just below the magnetic atom, 0 
is the angle formed by the N-TM-N trimer and the red arrow 
shows the position of the Cu atom used to define the surface. 

the TM introduces a local distortion on the CU 2 N lattice, 
shown in Fig. [2] the underneath Cu atom is pushed 
towards the bulk, and the N atoms are pulled out. These 
results are in line with previous worksi ^^i^^'^^ In order 
to characterize this structural distortion, we introduce 3 
distances: the TM- N distance {(In), the TM-Cu distance 
(dcu) and the TM-surface vertical displacement, z, that 
we take as the 2 component of the vector that joins the 
TM with the farthest Cu surface atom, marked with an 
arrow in Fig [21 In addition we also introduce the angle 
formed by the N-TM-N trimer (6). 

Our calculations, performed both for the small and 
large supercell (see Fig. [T|) show how these structural 
values are similar for different TM atoms, but with clear 
and systematic variations as a function of the number of 
d electrons. Whereas both the TM-N and the TM-Cu 
distances undergo minor variations across the TM series, 
the 9 angle has a much more marked change, going from 
structures where the TM is clearly a protrusion and the 
N atoms are weakly detached from surface, for small qd, 
to structures where TM is almost collinear with the N 
atoms, for Co and Ni. The tendency to form collinear 
N-TM-N structures is particularly clear in the case of Co 
and Ni chains (small cells). In the case of Co, the marked 
difference between the small and big cells suggests that 
there is a cooperative distortion in the case of Cobalt 
chains along the N direction,— also visible for Ni, and 
clearly absent in the case of the lighter TM, such as Fe, for 
which these chains have been studied experimentallyiii^ 

C. Electronic Properties 

We now discuss the electronic properties of the ad¬ 
sorbed TM on the Cu site of Cu2N/Cu(100). Our cal- 



FIG. 3. (Color online) Comparative study of structural prop¬ 
erties of TM atoms on CU 2 N. (a) Distance form the TM atom 
to the nearest neighbors N atoms d]v as a function of the TM 
atom, (b) Distance from the TM atom to the Cu atom lying 
just below the magnetic atom {dcu) ■ (c) Distance from the 
TM atom to the CU 2 N surface (z) as a function of the TM 
atom, (d) Angle formed by the N-TM-N trimer. Blue (red) 
line shows the result with the small (big) cell. 

culations for the pristine surface show that the Cu and 
N atoms in the CU 2 N layer have charge gcu = +0.2 and 
gN = —0.4. The adsorption of the TM atom results in a 
charge transfer mostly from the TM s orbitals to the N 
ligands, increasing their negative charge. The TM atoms 
lose practically all the 4s-electrons. The outermost elec¬ 
trons are thereby in the d shell. 

A naive interpretation of the picture that arises from 
the use of quantized spin models to describe these sys¬ 
tems would lead to conclude that the charge in the d shell 
of the TM is quantized. Our calculations show that this 
is not the case. A hint of this can be already seen by in¬ 
spection of the spin-resolved density of states projected 
over the d orbitals of the adsorbed TM atoms, shown in 
Fig. m obtained with QE, for the small cell. For in¬ 
stance, the occupancy of the majority spins (left-panel) 
is not 5 for Ni and Co. The presence of very broad peaks 
indicates strong hybridization of some of the d orbitals 
with the rest of the system, as we show below. 

In Fig. [5] we show the integrated density of states up 
to the Fermi energy, that gives the occupation of the d- 
shell for the different TM atoms, using both GCA and 
GGA+C, for several values of U. The results are plotted 
together with the charge of the isolated atom. In all 
instances we find that the d shell is more charged than 
in the free atom case. 

The difference between the computed charge and iso¬ 
lated charge, Ag, is shown in the inset for ?7 = 0 and 
U = 5eV. Expectedly, increasing U reduces Ag. Inter¬ 
estingly, Ag increases as we move away from half filling 
(the Mn atom). Varying U does not yield large changes 
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E-Ep [eV] E-Ep [eV] 

FIG. 4. (Color online) (a) DOS projected (PDOS) over d- 
orbitals of the different magnetic atoms, for majority spin 
(left panel) and minority spin (right panel). 



FIG. 5. (Color online) GGA+I/ calculation of the total charge 
of the 3d-levels for the different TM atoms and different values 
of U. Black dashed line show the result for the isolated atoms. 
Inset: deviation of the charge in the d levels with respect to 
the isolated atom 


in these results, except for the Mn atom, in line with re¬ 
sults obtained in reference |H for Co and Mn. We have 
also verified that Aq is stable with respect to changes in 
the size of the supercell. 


D. Magnetic Properties 

We now discuss the evolution of the magnetic mo¬ 
ments, fj,, of the series of 3d TM adsorbed atoms. In 
Fig. we show the magnetic moment of the free atoms, 
/ifree, as given by Hand’s. The largest free atom moment 
is n = 5fj,B, for the half-filled shell (Mn) and goes down as 
we move away from half filling. The upper panel of Fig [5] 
shows the GGA-hU calculation of the magnetic moment 
of the TM atoms. With the only exception of Cr, the 
magnetic moment for the adsorbed TM is always smaller 
than the free atom case. The deviations become par¬ 



Ti V Cr Mn Fe Co Ni 


• U=0.0 eV 
■ U=1.0eV 
*—• U=2.0 eV 
*—* U=3.0 eV 
X—K U=4.0 eV 
*—• U=5.0 eV 
•- -• Free Atom 


FIG. 6. (Color online) GGA-I-I/ calculation of magnetic mo¬ 
ment of the TM atoms (in units of ps), as a function of U. 
(a) Atomic (b) Unit cell. In both panels the black dashed line 
shows the expected magnetic moment for the free (isolated) 
atom (^free). 



FIG. 7. Orbital breakdown of the spin polarization, calculated 
using the small cell and f/ = 0. 


ticularly severe as we move away from half filling. For 
instance, the magnetic moment of adsorbed Ni is half 
the value of the free atom case. In contrast with the case 
of Ag, the value of the magnetic moment depends more 
strongly on the value of U, yet, the atomic limit is only 
reached in the case of Cr for U = 5eV. 

Whereas most of the spin is localized on the magnetic 
atoms, a substantial amount of the spin density is not 
located at the d-levels of the atom. Therefore, we also 
plot the cell magnetic moment and compared them with 
the free atom case (see Fig. [6 |d). It is apparent that, in 
the case of the cell moment, the deviations from the free 
case are very small for V, Cr, Mn and Fe. In the case of 
Ti and Ni, the deviation from the free case is still a factor 
of two. The overall trend is that, close to half-filling, the 
DFT results are closer to the free case. 

Fig. [7] shows the Spin polarization of the different d 
orbitals as a function of the TM atom. Whereas at half 
filling the magnetic moment has to be evenly distributed 
in the 5 d orbitals, away from half filling this is no longer 
the case. The orbital composition of the magnetization 
is interesting because it affects the magnetic anisotropy 
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Spin polarization of TM atoms 

Atom 

d-shell 

TM atom 

Cell 

Free atom 

Ti 

0.62 

0.71 

0.6 

1 

V 

1.32 

1.45 

1.41 

1.5 

Gr 

1.98 

2.10 

1.97 

2 

Mn 

2.25 

2.36 

2.47 

2.5 

Fe 

1.62 

1.75 

1.91 

2 

Go 

1.05 

1.1 

1.25 

1.5 

Ni 

0.46 

0.47 

0.57 

1 


TABLE I. DFT calculation of the spin polarization for {U = 5 
eV). 



FIG. 8. (Color online) GGA+t/ calculation of the magnetic 
moment of the nearest neighbor N atom as a function of the 
TM atom, for different values of U. 


and because only one of the d orbitals, the dl, couples to 
the s orbital of the last atom in STM tip. This results 
particularly interesting in the case of Ti, Co and Ni where 
we can appreciate that the d^ orbital has a very small spin 
polarization. These calculations were performed for the 
small cell, then we should expect interesting features in 
STM experiments when dealing with Co and Ni chains»ii 

The difference between the atomic and the cell mag¬ 
netization, summarized in Table HI implies that the sur¬ 
rounding atoms gain some magnetic moment as well. Our 
results, shown in Fig. [8] also show that their alignment 
with the N atoms can be both ferromagnetic (FM) or 
antiferromagnetic (AFM) depending on the TM. In par¬ 
ticular, the correlation is AFM below half-filling and FM 
above half-filling. At half filling the results depend on 
the value of U. Below we provide an explanation to this 
sign, based on a multi-orbital Anderson model. 

Finally, calculations for Co atoms were performed us¬ 
ing ELK for different values of U in order to check the 
magnetic properties obtained with QE. The results ob¬ 
taining with both codes were in very good agreement. 


III. CONNECTION WITH THE SPIN MODELS 

The DFT results of the previous section clearly show 
that the charge and spin of the d shells are not quan¬ 
tized. There is thus an apparent conflict between the 
DFT calculations and the use of spin Hamiltonian mod¬ 
els with quantized spinsfii^d^ In this section we address 
this important topic and provide a solution for this co¬ 
nundrum. First, we build a multi-orbital Anderson model 
starting from the DFT results. The Anderson model de¬ 
scribes the d orbitals of the TM atom, their hybridiza¬ 
tion to their neighbors, the crystal field splitting due to 
electrostatic interactions, intra-shell Coulomb repulsion, 
and the spin-orbit coupling. Most of these parameters 
are obtained from the DFT calculations, as we discuss 
below. Second, we solve the Anderson model exactly 
within a restricted multi-particle Hilbert space that in¬ 
cludes both and configurations, where d” 

stands for n electrons in the d shell of the TM and the 
p™ stands for m electrons in the p shells of the Nitro¬ 
gen first neighbors (without charge transfer there are 6 
electrons in each p-shell of the N atoms). Charge fluctu¬ 
ations were shown to be important in the case of Cobalt 
adatoms on MgO/ Ag<^ Including charge fluctuations in 
the Anderson model gives rise to a non-integer occupa¬ 
tion of the d shell but still preserves many-body states 
with a quantized spin S, that is identified with the spin 
of the quantum spin Hamiltonians. 


A. Maximally localized Wannier functions as 
atomic like basis set 

The derivation of an effective Anderson model start¬ 
ing from the DFT calculations requires a representation 
of the Khon-Sham Hamiltonian in a basis set that con¬ 
tains atomic-like d states localized around the transition 
metal atom. Our DFT calculations are performed with 
a plane-wave basis whereas the multiorbital Anderson 
Hamiltonian demands a local basis. Thus, to go from 
plane wave to a local basis, we represent the DFT Kohn- 
Sham Hamiltonian in the basis of maximally localized 
Wannier functions (MLWF))^“— computed using the 
code WannierQO, as described in reference Is^. 

The computation of the MLWF is implemented as fol¬ 
lows. First, we select a group of Bloch bands from a spin 
unpolarized^ calculation for a given TM/CU 2 N system. 
An energy window of 16eV around the Fermi energy is 
taken, and the band disentanglement procedure is per¬ 
formed. The selected Bloch bands are initially projected 
over the s, p and d orbitals of both the TM atom and the 
copper atoms and over the p and s orbitals of the nitro¬ 
gen atoms. The total number of states involved in this 
procedure is 313, corresponding to the 4 Nitrogen atoms 
(4 orbitals each), 32 Copper atoms and the TM atom (9 
orbitals each). 

An iterative procedure yields a total of 313 MLWF 
including 3 p-like and 5 d-like MLWFs localized around 
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FIG. 9. (Color online) Wannier orbitals for Fe and Co atoms 
at Cu 2 N. (a) py-orbital of a nearest neighbor N atom for Fe 
at CU 2 N close to the atom (isosurface 3). (b) p^-orbital of a 
nearest neighbor N atom for Fe at CU 2 N far from the atom 
(isosurface 1). (c) py-orbital of a nearest neighbor N atom for 
Co at CU 2 N close to the atom (isosurface 3). (d) p^-orbital of 
a nearest neighbor N atom for Co at CU 2 N far from the atom 
(isosurface 1). 


FIG. 10. (Color online) Wannier orbitals for Fe and Co atoms 
at CU 2 N. (a) dj. 2 _y 2 -orbital of a nearest neighbor N atom for 
Fe at CU 2 N close to the atom (isosurface 3). (b) d^ 2 _y 2 - 

orbital of a nearest neighbor N atom for Fe at CU 2 N far from 
the atom (isosurface 1). (c) d^2_y2-orhital of a nearest neigh¬ 
bor N atom for Co at CU 2 N close to the atom (isosurface 3). 
(d) d^ 2 _j, 2 -orbital of a nearest neighbor N atom for Co at 
CU 2 N far from the atom (isosurface 1). 


the Nitrogen atoms and the TM atom respectively. They 
are shown in Fig. [9] and Fig. [TOl where two different 
iso-surfaces of the MLWF with y and , for Fe 

and Co, are represented. In the case of iso-surfaces with 
larger value, corresponding to the wave-function close to 
atomic cores, these MLWF have the same symmetry as 
the real spherical harmonics with L = 1 and L = 2. In 
comparison, the iso-surfaces with small value do not show 
the symmetry of the Cartesian atomic orbitals. 

The Anderson model is build in a basis of single par¬ 
ticle states that involves the 5 d orbitals of the TM and 
the 6 p orbitals of the two first neighbor nitrogen atoms. 
The representation of the Khon-Sham Hamiltonian Hks 
in this basis can be written as: 

^CF + ffhyb= f (1) 

y ^pp j 

where we identify the crystal field Hamiltonian for the 
d orbitals dicF = Hdd and the pd hybridization Hamil¬ 
tonian with the off-diagonal blocks Hdp and Hpd- Im¬ 
portantly, Fig. [9] and Fig. [10] show how the MLWF for 
Cobalt are more extended than those of Fe, both for the 

—y^ and y orbitals. This accounts for the fact that the 
hybridization between the TM d orbitals and Nitrogen p 
orbitals is larger for Co than for Fe. In particular, the 
matrix element — y'^\HKs\y) is twice as large for Co 
than for Fe. 


B. Multi-orbital Anderson model 

We now introduce the multi-orbital Anderson model 
for the TM on the CU 2 N substrate describing the elec¬ 
trons in the 5 d orbitals hybridized with the p orbitals of 
the two adjacent N atoms. The Hamiltonian is the sum 
of 4 terms, the Kohn-Sham Hamiltonian of the previous 
section, plus the Coulomb repulsion of the electrons in 
the d shell, the Coulomb attraction of the d levels with 
the TM nucleus, and their spin-orbit coupling: 

H = Hks + -^Coui + Hz + Hso (2) 

In the following we label the five d-like MLWF of the 
TM with the index to, and 6 p-like MLWF of the two Ni¬ 
trogen atoms with the index n. The second quantization 
representation of the first term reads: 

Hks + Hz= + EdS m,m' )dlnad m'cr 

m, m' ,<7 

+ X! (^I^ksI«')pLp"'<7 + 

n, n' ,(7 

+ {{m\HKs\n)dl^^Pna + t^.c.) ( 3 ) 

m,n,cr 

The matrix Hdd = J2m,m 'describes the 
crystal fieldi^ The Ed energy scale accounts for the 
Coulomb interaction with the positive charge in the nu¬ 
cleus and has to be included to offset the excess in 
Coulomb repulsion in the configurations with an 

extra electron. The matrix Hpp = {n\HKs\n') describes 
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the single-particle p levels of the first neighbor Nitro¬ 
gen atoms and Hdp = describes their hy¬ 

bridization with the d levels, responsible of the charge 
fluctuation. Hdd, Hpp and Hpd are obtained from the 
DFT Hamiltonian using the wannierization procedure de¬ 
scribed above. 

The electron-electron Coulomb repulsion in the d shell 
reads: 


HCoul = 


E 


v; 


mm m m 




ma^m” <7'^^ 


( 4 ) 

For the evaluation of the Coulomb integrals Vmm"m'm"' 
we transform the angular part to a basis of eigenstates 
of = 2. For the radial part we take an effective radial 
hydrogen-like function ( with effective charge Z and a 
effective Bohr radius a^) to avoid the otherwise cumber¬ 
some numerical integration of the actual Wannier func¬ 
tions. In the basis of eigenstates of £ = 2, all the Coulomb 
integrals scale linearly with the value of Vbooo = 
Altought the numerical evaluation of U in terms of the 
parameters Z and is straight-forward, for the sake of 
generality, U will be considered a parameter, taken to sat¬ 
isfy the atomic Hund’s rule. For a given choice of f/, Ed 
is adjusted so that the average charge in the ground state 
of the many-body calculation described below equals the 
charge obtained in the DFT calculations. 

The spin-orbit term in the TM reads: 


Hso = Aso S\m'a')d\^„dm'G', (5) 

mm' ,<7(7' 


where \so is the atomic spin-orbit coupling of the d- 
electrons. This term is the only one that does not com¬ 
mute with the total spin operator and is the ultimate 
responsible of the lifting of the 2 S' -|- 1 degeneracy of the 
spin multiplets. 

In the following we shall show results for the case of 
Fe on CU 2 N. Given the small size of the single particle 
basis (5 d orbitals and 6 p orbitals) and the fact that we 
restrict the Hilbert space to the configurations d^p^^ and 
d”“l'ipii, with n = 6 for Fe, with a total of 1650 multi¬ 
electron states, the multi-orbital Anderson model can be 
solved by exact numerical diagonalization. 


C. Effective spin model 


In the seminal worb^ of Hirjibehedin et al., the spin 
excitations measured with STM- lETS were found to be 
described with the following spin Hamiltonian: 


H = D 


(eVs)‘ 


■E 




( 6 ) 


with 5 = 2 and Ci = (0,1,0) along the Nitrogen di¬ 
rection, and 6*2 and 63 are the off-plane and the hollow 
directions. The experimental results could be fitted with 


D = —1.55meV and E = 0.31 meV. Thus, the wave func¬ 
tion of the ground state and first excited state would be 
given by linear combinations of the rates of | 2 , ± 2 ), with a 
small mixing with the state | 2 , 0 ) in the case of the ground 
state. The height of the inelastic steps was found to be 
related to the matrix elements of the spin operators, giv¬ 
ing additional support to the notion that the quantized 
spin Hamiltonian ([S]) provides a quite good description 
of the spin excitations of iron on this surface. Recent 
experiments^^ with a detailed study of the lETS of sin¬ 
gle Fe/Cu 2 N as a function of the three components of 
the magnetic field show that the addition of extra terms 
in the Hamiltonian ([ 6 ]) yields an even better agreement 
with the experiment. Spin chains formed with Fe atoms 
in this system can also be modeled successfully with this 
Hamiltonian and the addition of interatomic Heisenberg 
couplingAltogether, these results support the notion 
that Fe can be described with a quantized anisotropic 
spin 5 = 2. 


D. Adiabatic continuity and spin conservation 
argument 


We now address the crucial question: given that ac¬ 
cording to DFT, both charge and spin of the d electrons 
in the transition metal are not quantized, what is the ori¬ 
gin of the quantized spin in the model Hamiltonian?. We 
now show that the quantized spin belongs to the many- 
body wave function that combines configurations with 
different charge states in the d shell. In the case of Fe, 
these would be configurations d^p^^ with = 2 and 
configurations with dJp^^ and = |, where is the 
spin of the d electrons. These many-body states yield 
non-integer charge and magnetic moment in the TM, in 
agreement with DFT, but they have a well defined to¬ 
tal quantized spin, in agreement with the spin quantized 
models. This is strictly true in the absence of spin-orbit 
interactions, and remains true when spin-orbit splittings 
are much smaller than the energy gap between different 
multiplets, which we show to be the case in a wide range 
of parameters. 

In order to understand the results of our numerical 
calculations, it is convenient to remind that atomic iron, 
with 6 d electrons, has a ground state with 5 = 2 and 
L = 2, and a total degeneracy of (2L -|- 1)(25 -I-1) = 25. 
The former ground state is captured by the multi-orbital 
Anderson model when switch off the Hdd crystal field, 
the Hpd hybridization, the SOC, considering configura¬ 
tions with 6 electrons. The crystal field Hdd quenches 
the orbital moment, so that in the absence of spin-orbit 
coupling, the ground state of the model has 5 = 2 and 
no orbital degeneracy. This multiplet is separated from 
the next higher energy multiplets, also with 5 = 2, by a 
gap of at least 300 meV, although this number depends 
on U . 

We now discuss the effect of mixing configurations with 
a different number of electrons in the d shell. On one 





FIG. 11. (Color online) Numerical diagonalization of the 
multi-orbital Anderson model for Fe atom at CU 2 N. (a) Aver¬ 
age charge of the d-levels of Fe atom for the Ground State as 
a function of Ed in atomic units {U = 5.5 eV and Aso = 50 
meV). (b) Sy {y is the N atoms axis) of the d electrons for 
the 5 lowest eigenvalues as a function of Ed {U ~ 5.5 eV and 
Aso = 0 ). (c) Excitation energies corresponding, for small 
Ed, to the 2 lowest multiplets with S' = 2, as a function of Ed 
{U = 5.5 eV and Aso = 50 meV) . 

hand, the relative weight of the configurations and 
depends on U, which is taken to he U = 5.5eV. 
On the other hand, i^pp, Hdd, Hpd, are obtained from 
DFT and Ed remains as the only adjustable parameter 
in the calculation. We start with a value of Ed so that 
the charge at Fe is g = 6 and ramp up Ed, so that q 
increases, as shown in Fig. [ffla). In Fig. me) we plot 
the evolution of the excitation energies as a function of 
Ed- For q = 6 the ground state has a multiplicity of 5. 
The next multiplet lies 300 meV above. As we ramp up 
q, these two multiplets remain well separated in energy, 
quite beyond the point where q = 6.4, the value obtained 
from DFT. We thus see that the ground state multiplet 
at g = 6.4 is adiabatically connected to the ground state 
multiplet at g = 6 . The effect of spin-orbit coupling, that 
we discuss in detail below, is to create a small splitting 
and to mix different states within the multiplet as Ed 
is varied, implying a change in the magnetic anisotropy 
tensor. 

The total spin is preserved as we ramp Ed- However, 


by changing the relative weight of d^p^"^ and dJp^^ config¬ 
urations, the magnetic moment in the d levels is expected 
to be reduced, moving from S{d) =2 towards S{d) =3/2 
. This is reflected in our calculations (see Fig. fTlT b)). 
taking Ago = 0 and a finite magnetic field By that lifts 
the 25" -I- 1 degeneracy. As Ed is increased, the expec¬ 
tation value of the operator describing the spin of the 
d electrons along the y axis, calculated with the 5 
lowest energy states, evolves from the eigenstates of Sy 
for S' = 2 to non-quantized values, in agreement with 
DFT results. 

Further increase of Ed yields that the mixing between 
d^p^'^ and d’^p^^ is so large that the splitting between the 
ground state S = 2 multiplet and the first excited mul¬ 
tiplet vanishes. When such regime is reached, the spin 
of the ground state changes, breaking the adiabatic con¬ 
nection with the state with S = 2 and quantized charge. 
In our calculations this happens for g cs 6 . 8 , larger than 
the DFT charge, g > 6.4. Thus, the model captures the 
crossover from the weak coupling limit, where the ground 
state adiabatically connected with the g = 6 ,5 = 2 , 
state, to the strong coupling limit in which the spin of the 
ground state multiplet changes and the adiabatic connec¬ 
tion is lost. 


E. TM-Nigrogen spin correlation 

The fact that the quantized spin corresponds to con¬ 
figurations with two charge states involving both d elec¬ 
trons in the TM and p electrons in the first neighbor 
Nitrogen atoms has implications on the spin correla¬ 
tion of the TM and N magnetic moments. Since both 
and have the same total spin, we have 

S'(d”) = St = S{d^~^^) ± i. The sign, and thereby the 
spin-correlation between the unpaired electron in the lig¬ 
and and the magnetic moment of the atom, depends on 
whether S{d'^) is larger or smaller than S'((i"+^). Thus, 
for Fe we have that S{d^) = 2 and S'(d^) = 3/2, so that 
the unpaired fermion must couple ferromagnetically with 
the S = ‘ij’l of the configuration, to keep 5 = 2. 
In contrast, for Cr we have 5(d^) = 2 mixing with 
S[d^) = 5/2 configurations, to bring up the charge, so 
that the unpaired fermion in the Nitrogen must couple 
antiferromagnetically with the atomic magnetic moment. 

This argument accounts for the trend obtained in our 
DFT calculations, shown in Fig. [ 8 l where the small mag¬ 
netization of the first-neighbor nitrogen atoms is antipar¬ 
allel to the TM magnetic moment for Ti, V, and Cr, and 
is ferromagnetic in the case of Fe, Co and Ni. In the case 
of Mn we obtain both signs, depending on U. From the 
argument of the previous paragraph we would expect a 
ferromagnetic coupling. Incidentally, the same argument 
can be applied to the conventional Anderson model with 
S = 1/2, predicting correctly the well knowit^ antiferro¬ 
magnetic interaction between the local moment and the 
adjacent electrons. 

























9 



FIG. 12. (Color online) Numerical diagonalization of the 
multi-orbital Anderson model for Fe atom at CU 2 N {U = 5.5 
eV and Xso = 50 meV). Eigenvectors projection over Sy for 
the ground state (upper panel) and the first excited state 
(lower panel) as a function of the charge in the d-shell. Blue 
lines show projection over |±2), red lines show projection over 
I ± 1) and green line shows projection over |0). Blue dashed 
line shows the projection over | ± 2) obtained from the spin 
model.— 


F. Spin-orbit coupling, magnetic anisotropy and 
symmetry of the wave functions 

The discussion above has ignored the role of spin- 
orbit coupling, even if the numerical results shown in 
Fig- HHa) and Fig. me) are obtained with a spin or¬ 
bit coupling Xso = SOmeV. These calculations show that 
the spin-orbit coupling splits the otherwise {2S -1-1) de¬ 
generate multiplets. However, the splittings between the 
different multiplets are still much larger than the energy 
gaps between them, so that a well defined total spin S can 
be atributed, even in the presence of spin-orbit coupling. 
The fine structure within the lowest energy multiplet can 
be described with an effective spin Hamiltonian, like the 
one in Eq. 

We now discuss the symmetry of the wave functions 
obtained from the exact diagonalization for two values of 
Ed, corresponding to having either exactly q = 6 elec¬ 
trons at the d shell or g = 6.4, the average charge ob¬ 
tained from DFT. For that matter, we use the fact that 
the eigenstates tpn of the multi-orbital Hamiltonian can 
be written as linear combinations of configuration states 
with well defined total Sy\ 

\^n) = 

Sy,l 

where 7 labels all the other quantum number necessary 
to characterize the basis set. In figures IT^ a.b). fT^ a.b) 



FIG. 13. (Golor online) Numerical diagonalization of the 
multi-orbital Anderson model for Fe atom at CU 2 N (U = 5.5 
eV and Xso = 50 meV). Eigenvectors projection over Sy for 
the second excited state (upper panel) and the third excited 
state (lower panel) as a function of the charge in the d-shell. 
Blue lines show projection over | ± 2), red lines show projec¬ 
tion over I ± 1) and green line shows projection over |0). Red 
dashed line shows the projection over | ± 1) obtained from the 
spin model.— 


and fT^ a) we plot the projection of the 5 lowest energy 
states of the multi-orbital Anderson model over the eigen¬ 
states of = 2,Sy). In the case of g = 6 it is appar¬ 
ent that the ground state wave function is dominated by 
Sy = 0 states, in disagreement with the experiment (see 
Fig. IT^ a)). The arrangement of the energy levels seems 
to indicate that the Nitrogen direction (y axis) is a hard 
axis in the problem (ei = (0,1,0) and Z? > 0 in Eq. ([G])). 
Given that the charge fluctuations are negligible in this 
limit, the main contribution to the magnetic anisotropy 
comes from the interplay between the crystal field term 
in the Hamiltonian, Hdd and the spin-orbit coupling. In¬ 
terestingly, when Ed is ramped so that g increases, the 
content of the wave functions evolves and for g = 6.4 
the wave functions (Eig. [T2l Fig. [13] and Fig. [Ha)) 
are in good agreement with those obtained from the spin 
model^ in which the Nitrogen direction is the easy axis 
in the problem (ei = (0,1,0) and Z? < 0 in Eq. ([ 6 ])). 
In particular, we note that the five lowest energy states 
of the Anderson model have wave-functions with strong 
overlap with those of the effective spin Hamiltonian<^ We 
thus see that the ligand field contribution, coming from 
the dp hybridization, changes qualitatively the magnetic 
anisotropy tensor. Interestingly, in the case of Fe/Cu 2 N 
we find that the inclusion of charge fluctuations is es¬ 
sential to capture the correct easy axis within the multi¬ 
orbital Anderson model. 

We finally analyze the low energy excitation spectrum, 
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FIG. 14. (Color online) Numerical diagonalization of the 
multi-orbital Anderson model for Fe atom at CU 2 N {U = 5.5 
eV and Aso = 50 meV). (a) Eigenvectors projection over Sy 
for the fourth excited state as a function of the charge in the 
d-shell. Blue lines show projection over | ± 2), red lines show 
projection over | ± 1) and green line shows projection over 
|0). Green dashed line shows the projection over |0) obtained 
from the spin model.— (b) Low energy spectrum as a function 
of the charge in the d-shell. 

En — Eq, where Eq is the energy of the ground state. In 
the range of Ed considered, such that q moves from the 
nominal value g = 6 to the DFT value q = 6.4, the 
five energy levels of the lowest energy multiplet are al¬ 
ways split. This is expected in the case of a integer spin 
{S = 2) described with Hamiltonian Eq. Interest¬ 
ingly, the low energy splittings increase as q is increased 
towards the DFT value, as shown in Fig. [T4l b). This 
is related to the fact that the energy gap between the 
first and second five-fold degenerate multiplet decreases, 
as shown in Fig. fTlT c). This behavior can be understood 
in terms of degenerate perturbation theory, where spin- 
orbit coupling only splits the states in the lowest energy 
multiplet through virtual transitions to the higher state 
multiplets. We note that the spin excitation energies ob¬ 
tained from our calculation are a 25 percent smaller than 
those observed in the experiment. So, whereas the model 
captures the right symmetry, it can only give a rough de¬ 
scription of the excitation energies, which is probably due 
to the approximations in the model, such as the restric¬ 
tions taken in the many-body Hilbert space. 

IV. CONCLUSIONS 

In this work we have undertaken a systematic study 
of the electronic properties of the 3d transition metals 
on the CU 2 N surface. We systematically find that the 


charge and spin of the d electrons are not quantized, 
and are thereby different from the one in isolated atoms, 
which is expected given the conducting nature of the sub¬ 
strate. We have then addressed the issue of how to recon¬ 
cile these results with the fact that quantized spin mod¬ 
els account for the spin excitations of Mn, Fe and Co 
ad-atoms. For that matter we propose a multi-orbital 
Anderson model in which many-body states that mixes 
configurations with two charge states in the d shell are 
considered. Importantly, even if the charge is not well de¬ 
fined in the d shell, these multi-electron wave functions 
have a well defined total spin S. We find that the states 
with a quantized charge in the d shell are adiabatically 
connected with the actual many-body states that mix 
configurations d^p^"^ and in the sense that both 

have the same total spin S and there is no mixing with 
higher energy multiplets as the addition energy is varied 
numerically. 

We thus conclude that quantized spin S of the model 
actually refers to the spin of these many-body states that 
include both the d electrons and the ligand electrons. It 
is thus fair to say that the magnetism in this system 
is not strictly atomic, which connects with previous re¬ 
sults in the case of magnetic atoms on metallic surfaces)^ 
for which sophisticated theoretical treatments have been 
proposed Our picture provides a natural explana¬ 

tion to the sign of the spin correlation between the TM 
and the nitrogen atoms obtained in the DFT calcula¬ 
tions and reconciles the use of quantized spin Hamilto¬ 
nians with the results of DFT calculations. Our Ander¬ 
son model calculations for Fe/Cu 2 N also indicate that a 
charge fluctuations in the Fe d shell, due to the hybridiza¬ 
tion to the ligands, are essential to capture both the sym¬ 
metry and the magnitude of the magnetic anisotropy ob¬ 
served in the experiment. 

Finally, we expect that our analysis should also be ap¬ 
plicable to other systems with magnetic adatoms and 
molecules deposited in conducting surfaces whose spin 
excitations can be described in terms of quantized spin 
model o^^'^^~ — 
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